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Abstract 

In the framework of a chiral constituent quark model, considering the contributions of tt anni- 
hilation and one-gluon annihilation, the proton-antiproton S-wave elastic scattering cross section 
experimental data can be reproduced by adjusting properly one-gluon annihilation coupling con- 
stant. Meanwhile, using the fixed model parameter, we do a dynamical calculation for all possible 
5-wave nucleon-antinucleon states, the results show that, there is no S-wave bound state as indi- 
cated by a strong enhancement at threshold of pp in J/ip and B decays. 

PACS numbers: 12.39.Jh, 14.20.Pt, 13.75.Cs 
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I. INTRODUCTION 



Theoretical studies of baryon-antibaryon bound states date back to the proposal of Fermi 
and Yang [l| to make the pion with a nucleon-antinucleon pair. The traditional NN inter- 
action studies, such as boundary condition model [2J, optical model [3| and coupled- channel 
models ^j], emphasized on the handling of the short-range part of NN interaction, and for 
the long range part they are quite similar. The extensive and excellent reviews are given 
in Ref . |5| . Although the agreement with scattering experimental data is obtained, from the 
quantum chromodynamics (QCD) point of view it is hard to image that a hadronic pic- 
ture can be applied to such a short-range where hadrons are "overlapped" and the internal 
structure of hadrons must be considered. 

A possible way out is to start from QCD, the fundamental strong interaction theory. 
Namely we should start with quark-gluon degree of freedom, rather than the meson-baryon 
picture. QCD has been already proved to be the right theory at high energy. At low energy, 
because of the non-perturbative nature of QCD, one has to rely on effective theories and/or 
QCD-inspired models to get some insight into the phenomena of the hadronic world. The 
constituent quark model [6, Q, 9] is one of them. It has been successful in describing 
hadron spectrum, the baryon-baryon interactions and the bound state of two baryons, the 
deuteron. Therefore, extending the constituent quark model to NN study is an interesting 
practice. 

The early studies in the traditional meson exchange framework found that, if neglecting 
annihilation channels, many bound states might be formed, while annihilation effects were 
included, the binding force decreased and some bound states were washed out [10]. Therefore, 
how to take into account the effect of annihilation in an unified framework is important. We 
attempt to include the contributions of the annihilations of a qq pair into meson or into gluon 
in an unified manner, besides including tt, a and gluon exchange in the constituent quark 
model. In order to keep the model well-describing baryon spectrum and NN scattering 
data, all of the model parameters have been fixed as much possible as those fixed in NN 
interaction and baryon spectrum. 

Perturbative QCD calculation showed that the gluon running coupling constant decreas- 



ing with the increasing of momentum transfer. Here we vary gluon annihilation coup! 
constant a' s to see if the proton-antiproton 5*-wave elastic cross section experimental data 
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can be well reproduced. For comparison, the case of ignoring the contributions of tt annihi- 
lation is also computed. 

The BES collaboration in the radiative decay JAfr — > 'jpp observed a sharp enhancement 
at threshold in the pp invariant mass spectrum [12]. They tried to fit the enhancement 
by means of a S-wave Breit-Wigner resonance, and obtained the resultant mass peak be- 
low threshold. Belle also reported they observed an enhancement in the pp invariant mass 
distribution near the threshold in the decays B + — > K + pp and B° — > D°pp 13J. Many inter- 
pretations Q, Q, 16, 17] on the observation were suggested. Here we apply the constituent 



quark model constrained by baryon spectrum and NN interaction and in addition consid- 
ering the contribution of annihilation fixed by proton-antiproton 5*-wave elastic scattering 
cross section, to do a dynamical calculation for all possible S'-wave nucleon-antinucleon 
system to study if there is pp S'-wave bound state. 

The paper is organized as follows. Sec. II explains the model Hamiltonian, its parameters 
and the calculation method. Sec. Ill is the results and discussions. 



II. HAMILTONIAN, MODEL PARAMETERS AND CALCULATION METHOD 

We take the chiral quark model used in the study of multi-quark system, which essentially 
is an effective theory on exchanges of Goldstone boson, scalar meson a, as well as gluon 
between quarks, and extend it to include the antiquarks and the annihilation interaction, to 
study the nucleon-antinucleon system. 

As the first step, here only .S'-wave states of the nucleon-antinucleon pair are considered, 
i.e., the total orbital angular momentum L = 0, and we have J = S (the total angular 
momentum comes from quark spin only). 

We start from a Hamiltonian which was used by Salamanca group for NN interaction 8] . 
It is one of the chiral quark model with only tt and a mesons. The quark-meson coupling 
constant a c h is fixed by gjf Nn /ATT 

9 2 ch ^ ,2 9 2 nn, m l „ 9ch m l 



4tt V 4tt M 2 n ' Ait Am 2 u ' 

The quark mass m q and a meson mass m a are taken to be 313 MeV and 3.421 fm _1 , 
respectively. Quark-gluon coupling constant a s is determined by the A — N mass difference. 
The confinement strength, a c is obtained from nucleon stability condition. Parameters and 
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the calculated deuteron properties are listed in Table I. 

Table I. Model Parameters, deuteron properties. 8oge, &ope (MeV) are the contributions 
of gluon, it exchanges to A — N mass difference, respectively. 



fe(fm) a s a ch A(fm 1 ) 


a c (MeV fm~ 2 ) 


Fo(MeV) 


6oge 5ope B D (MeV) vTr 2 )(fin) 


Pd 


0.518 0.485 0.027 4.2 


46.938 


-487.29 


145.6 148 2.0 1.96 


4.86 



To extend model from NN systems to NN systems, we have to take into account the 
annihilation contributions in addition to the scattering ones: the relevant scattering and 
annihilation Feynman diagrams are shown in Fig.l and Fig.2. 
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FIG. 1: The exchange diagram between quark and antiquark 
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FIG. 2: The annihilation diagram between quark and antiquark 

Taking n exchange as an example (in orbit and spin space), the T-matrix of tt exchange 
diagram between quark and antiquark can be written as 

— F 2 ((f) 

u(Pv s'lhVPi, s i)^- t T9ch^(P2, s 2 h 5 v(p' 2 , s' 2 ), (1) 

q z + 7774 



here, a form factor F(cp) is assumed to be (A 2 /(A 2 + q 2 )) 1 / 2 . p lf Si,p[, s[ are the four- vector 
momenta and spin z-projections of initial quark and final quark, respectively. p 2 , s 2 ,p 2 , s 2 
are the four-vector momenta and spin ^-projections of initial antiquark and final antiquark, 
respectively, u and v are assumed to be free Dirac spinors for quarks and antiquarks, 
respectively, is the ir mass. 

Similarly, T-matrix of it annihilation diagram can be written as 

-1 



-M(p' 1 ,s' 1 )7 5 f(p2,4 



2 q 2 + ml 



(2) 



using Fierz identities for Dirac matrices, a four-fermion matrix can be expressed as a linear 
superposition of other matrices with a changed sequence of spinors 18]. 
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lOibXcO'd) = J2(aO k d){cO k b) 



k=l 



Then we take the non-relativistic limit and transform the potential to space-time repre- 
sentation, we have 
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(4) 
(5) 



and Vg d a are the effective potential from 7r exchanges and annihilation between quark 



gg 



and antiquark, respectively. V q q are the effective potential from one gluon annihilation. 
Because here N and N are both color singlet and there is no quark exchange between 
and N, one gluon exchange between N and N does not contribution at all. 

Since we are taking into account the contributions of the lowest order to the S'-wave at this 
step, under static approximation the contribution from a meson annihilation in the present 
case vanishes. The a exchange potential between quark and antiquark in the non-relativistic 
limit in coordinate space can been written as 



gg v ; 



9ch 



A 2 



mi 



-Ar 



(6) 



An A 2 

The detailed derivation of these effective potentials between quark and antiquark can be 



found in Ref. 



19J. 
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Therefore, the Hamiltonian of NN system is 

H pp = EK + ^~)~ Tcm + E [V conf (n 3 ) + V e (r l3 ) + V^( nj )] 

1=1 % 2>J=1 

V^nj) = V£ >qq ( nj ) + V£ qq ( rij ) + V^{ Tij ) (7) 
V&rv) = V^{r l3 ) + V°({r l3 ) + V^{r l3 ) 
V conf (r i3 ) = -a c A- • XjTij + V . 

We should note that, in order to keep the model well-describing baryon spectrum and NN 
scattering data, all of the model parameters related to NN system are unchanged, only one 
parameter a' s connected with annihilation has been left to be adjusted. 

We use Kohn-Hulthen-Kato variational method for bound and scattering problems in an 
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unified framework 
(i) for bound problem 

Following the cluster model approach, the resonating group method (RGM) wave function 
is written as 

here R is the relative coordinate between the clusters of B\ and B2. A is the antisym- 
metrization operator but in fact there is no need for this antisymmetrization for NN system 
because q and q can be treated as different particles. (f>B 1 and (f>B 2 are the internal wave 
functions of two quark clusters, x{R) is relative motion wave function. 
The relative motion wave function is expanded into partial waves 

N 

X (R) = £ x L {R)Y LM {R) = E E cf xf (R, SjY^iR) (8) 

L L i=l 

with 

= M ^ )3/4expi ~^ {R2 + S ^^ RS * )YLM{R) (9) 

where ii is the modified spherical Bessel function. In this paper, only S wave (L = 0) is 
taken into account. 

Adding the center of mass motion, the wave function of six quarks can be written as the 
production of the single-particle orbital wave function with different reference centers, i.e., 

n „ 3 6 

L k i=l a=l 13=4 



here ip a (Si) and i^^—Si) are the single-particle orbital wave function with different reference 
centers, is the coupled channels index. 

Via the variation with respect to the relative motion wave function the RGM 

equation 

J H(R, R')x(R')dR' = E J N(R, R')x(R')dR' (11) 
becomes an algebraic eigenvalue equation 

Y.C h uH-f k = E^C^N^. (12) 

j,L' k j 

where N { J and k are the wave function overlaps and Hamiltonian matrix elements, 

respectively. 

(ii) for scattering problem 

The wave function of the relative motion is expanded by 

n 

X (R) = £ X \R)Y LM {R) = £ E cf SdY LM {R) (13) 

L L i=0 

here 

m,s,) = { () fl< ^ (u) 

[ {hl\kR) - Si hY\kR))kR R> R c 

with 

where il is the modified spherical Bessel function, is the L-th spherical Hankel function, 
k = ^f(2fiE re i). 

The constants «j and Sj are determined by the condition that the relative motion wave 
function for R < R c and R > R c smoothly connect at R = R c . 

Using the normalization condition cq = 1 — X^=i c « an d varying with respect to parameters 
Ci (i = 1, ...,n), we have 

([tij - M ]L L M t ) = 0,i = 1, ...,n 

with a symbol {fh^g) by 



/*oo /*oo 

(fL L g)= / f(R)L L (R, R')g(R')dRdR' 
Jo Jo 
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L L (R,R') = J [<j>$(U)<f>ttZB)uf(R)Y£ M (R)](H - E)A[MU)M^)uf(R)Y LM (R)], 
then we have the n linear equations for the q's, 



where 



with 



5^ hijCj = Mi (16) 



Ly = Kij - K i0 - K oj + K 00 , M, = K 00 - K i0 (17) 



Kij = (uiL L Uj) (18) 



The final phase shift is 



n 



S st = S t + iaJ2 K 0i c i ( 19 ) 

8=0 

here a = fj,/k, // is the reduced mass. The difference s t — s st is a good measure to check 
the accuracy of the calculation. In order to calculation (term(18)) as analytically as 
possible, there is a very useful skill mentioned in Ref. 



III. RESULTS AND DISCUSSIONS 



There are four possible states for an S-wave nucleon-antinucleon systems with different 
isospin and total angular momentum J respectively. They are IJ PC = 11 , 10 _+ , 01 , 
00 h . Since running coupling constant would change with momentum transfer, here we 
adjust gluon annihilation coupling constant a' s from being smaller than exchange coupling 
constant a s to larger than a s . 

In experiment one measures the spin-averaged scattering amplitudes, and the elastic 
scattering amplitude is the sum of the isospin 1 = and 1 = 1 amplitudes with equal 
weights. According these, we calculated the total proton-antiproton elastic cross sections, 
and compare it with experimental data. 

We calculated pp S-wave elastic cross section including 7r annihilation and gluon annihi- 
lation with different gluon annihilation coupling constant a' s , the results are shown in Fig.3. 
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We found that, if we let a' s = (i.e., we did not take into account the effect of gluon anni- 
hilation) the cross section would be very larger than experimental data. It implies that, if 
excluding the gluon annihilation process, NN system is more attractive than NN system 



and many bound states might be formed, which is consistent with the results in Ref.|10l. 122] . 
When we add the contribution of gluon annihilation, with a gradual increasing a' s , the cross 
section will decrease quickly, especially in the low energy region. And the larger a' s is, the 
smoother the change of cross section with the scattering energy is. If we choose a' s equal 
to a s , the cross section will be below the experimental data. When we choose a' s is about 
the one third of gluon exchange coupling constant a s (a' s = a s /3), the cross section became 
close to experimental data, the difference between theoretical cross section and experimental 
one would not more than 5 mb. 
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Fig. 3 proton-antiproton elastic cross sections, including ir annihilation and one-gluon 
annihilation with different coupling constants. The full squares show experimental data of 
total elastic cross section. The open squares are that of S'-wave component [11]. 
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Fig. 4 Same as Fig. 3 but it annihilation excluded. 

Fig. 3 proton- antiproton elastic cross sections, including ir annihilation and one-gluon 
annihilation with different coupling constants. The full squares show experimental data of 
total elastic cross section. The open squares are that of S- wave component Since 
annihilation process would happen at short range, and 7r term might play an important 
role only at long range, so we do a calculation by excluding the n annihilation term. Fig. 4 
shows the pp cross section with the different a' s in the case of excluding the contribution of 
7r annihilation. We found that, larger a' s corresponds to smaller cross section. In order to 
reproduce pp .S'-wave elastic cross section, the effective gluon annihilation coupling constant 
a' s should increase to compensate for the absence of ir annihilation. The proper value is 
a' s =a s /0.9, which is very close to gluon exchange coupling constant. 

The above results implies that both gluon annihilation and ir annihilation provide effective 
repulsion, which would decrease the S'-wave cross section. By adjusting the strength of 
annihilation term properly, proton-antiproton S'-wave cross section experimental data can 
be reproduced no matter the 7r annihilation is included or not. 

Taking a' s = a s /3 and a' s = a s /0.9 as examples, we give the contributions of different 
spin and isospin components to the total cross section in Fig. 5 and Fig. 6 for including and 
excluding tt annihilation, respectively. The results show that, the contribution of I J = 00 
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channel is always very small no matter the tt annihilation is included or not; the I J = 01 
channel contribution to cross section is larger than I J = 10 channel when tt annihilation 
is included; while the former contribution is smaller than the latter if the tt annihilation 
is excluded; and with the increasing of scattering energy the difference between these two 
channels contributions are both decreasing in these two cases; For the case of including tt 
annihilation, the total cross section is close to the I J = 11 channel ones when scattering 
energy E cm is below 15 MeV and it will be close to I J = 10 channel ones when E cm is above 
15 MeV. If excluding the contribution of tt annihilation, the total cross section is very close 
to IJ = 11 channel in the whole energy range. 
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Fig. 5 S'-wave elastic cross sections with a' s = a s /3.0 in the case of including tt annihilation, 

for channel I J = 11, 10, 01, 00, respectively. 
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Fig.6 Same as Fig.5 for the case of excluding 7r annihilation with a' s = a s /0.9. 



The effective potentials of I J = 11, 10,01,00 channels are given in Fig.7 and Fig. 8 cor- 
responding to different annihilation coupling constants. From these figures we found that 
I J = 00 channel has very small intermediate attraction no matter what parameters are cho- 
sen, which is consistent with the results of cross section in Fig. 5-6. The other three channels 
all have attraction to some extent. Excluding the contribution of tt annihilation there will 
be more attraction left in the other three channels, and the minimum values of potential are 
all at the separation of 1.0 fm. 

The BES collaboration in the radiative decay J/ip — > ^ypp observed a sharp enhancement 
at threshold in the pp invariant mass spectrum[l2j, and they obtained the mass peak below 
the threshold of pp. The above effective potentials give qualitative information only, in order 
to see whether there is a pp bound state, we do a dynamical calculation. Taking into account 
the contribution of annihilation and using the model parameters (in Fig. 5-6) determined by 
proton-antiproton S-wave elastic scattering cross section, we solve RGM equation for all 
possible S-wave nucleon-antinucleon states, to see if there is pp S-wave bound state. Here, 
we use 10 basis wave functions to expand the wave function of the relative motion, the 
boundary point is at 5.8 fm. We find that, although the four possible channels with different 
isospin and spin quantum numbers all have intermediate range attraction ( from several Mev 
to about 50 MeV), there is no bound state in dynamical calculation. Moreover, there is no 
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resonant state below the system mass of 3.9 GeV in our calculation. That is to say, if we 
determine the model parameter of quark- ant iquark annihilation by fitting proton-antiproton 
S-wave elastic cross section experimental data, our model does not support a tight bound 
state claimed by BES and Belle experimental groups. 




-50 - 

.100 I 1 1 1 1 1 1 1 1 1 1 1 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 



S (fm) 

Fig. 7 Proton-antiproton effective potential with a' s = a s /3.0, including n annihilation, for 

channel I J — 11, 10, 01, 00, respectively. 

To sum up, in the framework of chiral quark model, adding the contribution of one 
gluon and n annihilation, we can reproduce proton-antiproton S-wave elastic cross section 
experimental data by adjusting the coupling constant of gluon annihilation term. Using 
the model parameters determined by pp scattering cross section, our dynamical calculation 
for NN system with IJ — 11, 10, 01, 00, quantum numbers does not find an S-wave bound 
state. 
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Fig. 8 Same as Fig. 7 for excluding ir annihilation with a' s = a s /0.9. 

Obviously our conclusion of no pp bound state is based on the assumption that the chiral 
constituent quark model is suitable for NN system. Our conclusion is also based on the 
assumption that the multi n channel coupling effect, which is possible within the chiral 
constituent quark model, can be neglected. In addition all of the hidden color channels 
coupling effects have been omitted. All of these effects should be studied further. 
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